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We calculate for the first time the electric dipole moment (EDM) of the ®Li nucleus within the 
a + p + n three-body cluster model using the Gaussian expansion method, assuming the one meson 
exchange P, CP-odd nuclear forces. It is found that the EDM of the ®Li is 2 times more sensitive on 
the isovector pion exchange P, CP-odd nuclear force than the deuteron EDM, due to the CP-odd 
interaction between the nucleons and the a cluster. The ®Be EDM is also calculated in the same 
framework as an a -|- a -|- n three-body system. We also test the ab initio calculation of the EDM 
of the deuteron, and ®He nuclei using the realistic Argonne ul8 nuclear force. In the ab initio 
calculations, good agreements with previous studies are obtained. We finally discuss the prospects 
for the new physics beyond the standard model. 
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I. INTRODUCTION 

It is commonly believed that the current baryon num¬ 
ber asymmetry was created at the early stage of our Uni¬ 
verse. It is however known that the standard model of 
particle physics is not able to realize it, and CP violation 
beyond the standard model (BSM) is therefore required. 
One promising experimental approach to search for the 
CP violation is the measurement of the electric dipole 
moment (EDM) [IH^. The EDM is an observable very 
sensitive to it, and has many advantages, such as the ac¬ 
curate measurability in experiments, low cost in the ex- 
perimental preparation, very small standard model con¬ 
tribution |10l |. etc, and it was so far measured in various 
systems, such as the neutron [n[, atoms [13, molecules 
[l^ . and muon [l^ . 

Among them, the nuclear EDM is of particular interest 
im. Recently, the measurement of the EDM of light 
nuclei using storage rings is planned at the Brookhaven 
National Laboratory, and the prospective sensitivity is 
O(10“^®)e cm [T 3 . The nuclear EDM also has its own 
advantages. First, the nuclear EDM does not suffer 
from Schiff’s screening [3 , since we have no electrons 
to screen the charged nucleus. Also, the nuclear sys¬ 
tem may enhance the nucleon level CP violation due to 
the many-body effect [l^l- These arguments make the 
nuclear EDM to be a very attractive probe of strong sec¬ 
tor CP violation, and theoretical investigations of the 
deuteron [l9l - [^ and three-nucleon systems psl-f^ were 
extensively done so far. 

As a natural choice of the next target of the discussion, 
we have the ®Li nucleus, since it is the lightest stable 
system with non-zero angular momentum which can be 
found after the three-nucleon systems. The ®Li nucleus is 
known to have a cluster structure [13 , so we may expect 
some enhancement of the nucleon level CP violation, due 
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to the derivative interaction of the CP-odd Hamiltonian. 
In this sense, the analysis of the ®Li EDM may also be 
interesting in the point-of-view of the study of the nuclear 
structure and the cluster dynamics M- 

In this work, we therefore investigate the EDM of the 
®Li nucleus in the cluster approximation, as the a—p — n 
three-body system. We also calculate the EDM of the 
®Be nucleus, which is calculable in the same framework 
[13, and a test of ab initio calculations of the deuteron, 
the ^He, and nuclei using the Argonne ul8 interaction 
[13 . To solve the many-body Schrodi^er equation, we 
use the Gaussian expansion method |34| . This method 
was applied in a wide number of subjects, extending from 
the particle to atomic physics (35l - [38l| . and it is also ex¬ 
pected to give accurate results in the study of the EDM 
of few-nucleon systems. 

This paper is organized as follows. We first give the 
definition of the EDM. We then briefly review our cal- 
culational method. The result of the EDM of the ®Li 
nucleus is then given, together with the EDMs of ®Be, 
^He, ^H, and nuclei. We finally see the prospect for 
the determination of the new physics BSM and summa¬ 
rize our discussion. 


II. THE NUCLEAR ELECTRIC DIPOLE 
MOMENT 

To induce the nuclear EDM, the existence of the P, 
CP-odd nucleon level processes is required. The nuclear 
EDM has two leading sources: (1) the intrinsic EDM of 
the constituent nucleons, and (2) the P, CP-odd N — N 
interactions (nuclear force) to polarize the whole nucleus. 
In this work, we neglect the exchan ge c urrent, since its 
contribution is expected to be small [20j. 

We hrst give the contribution of the intrinsic EDM 
of the constituent nucleons to the nuclear EDM. As the 
nucleon EDM is proportional to the nucleon spin, the 
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effect is simply given by 

i 

— {^p)a dp -j- {o'n)A dm (1) 

where | A) is the polarized nuclear wave function (A = 
^H, ^He, ®Li, ®Be). Here dp and d„ are the pro¬ 
ton and neutron EDMs, respectively. These are given 
parameters which depend on the QCD and elementary 
level physics. The coefficients {(Jp)A and (cr„)^ are just 
the spin matrix elements of the nucleus, and depend only 
on the nuclear structure. The contribution of the single 
nucleon EDM may be enhanced if the nucleon is relativis¬ 
tic inside the nucleus, as for the atomic systems [^.13.13^. 
As the nucleons are nonrelativistic in light nuclei, the lin¬ 
ear coefficients of the systems in question will not receive 
sizable enhancement. 

The effect of the nuclear polarization generated by the 
CP-odd nuclear force may, in contrast, be enhanced even 
in light nuclear systems. The polarization contribution 
of the P, CP-odd nuclear force to the nuclear EDM is 
given by 


A 

dT^ = ^|(i|(l + r;)7^,,|i), (2) 

i=l 


isovector, and isotensor structures. The Yukawa function 
is given by yx{r) = . As the CP-odd effect is 

small, the polarization contribution to the nuclear EDM 
is given by the linear term 

■ir’=E“«GS'. (4) 


The linear coefficients cS^x depend only on the nuclear 
structure, and are the main targets of this discussion. 


III. METHODOLOGY AND MODEL SETUP 


We study the structure of the ®Li and ®Be nuclei within 
the framework of the a + n+p and a -I- a -I- n three-body 
cluster models. Here we assume the a cluster to be an 
inert core. Also, we take three kinds of Jacobian coordi¬ 
nates for those systems (For example, see Fig. 3 of Ref. 
[4l[ in the case of the ®Li). Regarding the calculation of 
three-nucleon systems such as the ^H and ^He nuclei, we 
take the Jacobian coordinates shown in Fig. 1 of Ref. 


The Schrodinger equation is given by 


( 5 ) 


where | A) is the polarized (in the 2 -axis) nuclear wave 
function, and r/ is the isospin Pauli matrix. TZiz is (the 
2 -component of) the position of the constituent nucleon 
in the nuclear center of mass frame. This permanent po¬ 
larization effect is realized through the mixing of opposite 
parity states. 

We now show the detail of the CP-odd nuclear force, 
necessary to generate the polarization. In this study, 
the CP-odd nuclear force is given by the standard one, 
based on one-meson exchange potential [j^. The CP- 
odd hamiltonian is given by 

HfT = ^ 

+ri y^{T) - G® ypir)] 

+ ir;^_-V[GW J.(r)-G«%(r) 

+ irI<T+-V[GW3;.(7-) + G«%(r) 

+GWj^p(r)-G«3;.(r)] 

+ - D • T 2 ) • V [G(, 2 ) y^{r) 

-Gf 3^p(r)]|,(3) 


where G^^ = gxNNd^xNN coupling constant of the 

CP-odd nuclear force with the exchanged mesons X = 
Tr,r],p,u!. The index i = 0,1,2 denotes the isoscalar, 


where 


i/ = T-f^t4b + Ppaulh (6) 

a,b 

with the kinetic energy operator T, and I4b is the interac¬ 
tion between constituent particles a and b [including the 
CP-odd hamiltonian of Eq. ([S])]. The OCM projection 
operator VpauU is given below for the ®Li and ®Be. The 
total wave functions for the ®Li and the three-nucleon 
systems are described in Eqs. (3.2) and (66) in Refs. [4l| 
and [ 3 ^ , respectively. The total wave function of the ®Be 
nucleus is described as 


(®Be) 

c=l nl,NL IK sS 


X (u)^/’]vi (R-c)] A Xi (^1)] 


/,(=) 


J J 


( 7 ) 


Here the operator Sa stands for the symmetrization be¬ 
tween the two a clusters. The spin function of the nu¬ 
cleon is denoted by x 1 . 

Followin g th e Gaussian Expansion Method (GEM) 
[ 13 ) US HE El] , we take the functional forms of 4’nim{'r), 
ilNLM{^) as 


(/>nim(i’) = r'e Ijm (r) , 

iPnlmCR) = , (8) 
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where the Gaussian range parameters were chosen ac¬ 
cording to the geometric progression: 

r„ = (n = 1 - n^ax) , 

Rn = RiA^-^ (A^=l-iV„,ax). (9) 

The angular momentum space l,L,K < 2 was found to 
be sufficient to obtain good convergence of the calculated 
results. 

The Pauli principle in the N — a and a — a systems is 
taken into account in the orthogonality condition model 
(OCM) [i^. The OCM projection operator Ppauii is rep¬ 
resented by 

14>auli = lim VA| 0 /(r„x))(</'/(r’ax)|, ( 10 ) 

X—¥oo 

f 

which rules out the amplitude of the Pauli-forbidden a—a 
and a — N relative states (j>f{i'ax) from the three-body 
total wave function [i^. The forbidden states are / = 
Os, Is, Od for X = a and / = Os for x = TV, respectively. 
The Gaussian range parameter b of the single-particle Os 
orbit in the a cluster (Os)'’’ is taken to be 6 = 1.358 fm 
to reproduce the size of the a cluster. We employ the 
a — N and a — a interactions so as to reproduce the 
scattering phase shift of the a — N and a — a systems at 
low energy Also we use the Argonne nl 8 nuclear 

force for the three-nucleon systems, and the Argonne v8' 
interaction for the TV — TV subsystem of the ®Li nucleus 

M- 

For the CP-odd nuclear force, we have used the folding 
Va-N{r) = J d^r'Vf>f{r-r')pa{r'), (11) 

where Vff is the radial function of the CP-odd nuclear 

force and Pa{r) = ^ 3 ^ 3 / 2 e~'^ is the nucleon number 
density of the a cluster normalized to 4. The nucleon 
number density was approximated by a single gaussian 
with the spread b = 1.358 fm. The folding of the CP- 
odd potential is shown in Fig. [1] It is important to 
note that the folding cancels the CP-odd nuclear forces 
for which the spin or isospin Pauli matrix acts on both 
interacting nucleons simultaneously. Therefore, only the 
isovector pion exchange, isoscalar rj exchange, isoscalar 
and isovector cj exchange CP-odd nuclear forces survive 
for the a — N system. For the a — a system, all CP-odd 
nuclear forces cancel. 

We must also note that the folding may overestimate 
the effect of heavy meson exchange r],p,u!, since the a 
cluster approximation partially averages the microscopic 
level physics. At the microscopic level, the nuclear force 
has a repulsive core which hinders the contribution of 
the short range CP-odd potential with heavy meson ex¬ 
change. This effect may break away when the folding 
averages the potential at the a cluster scale. We will re¬ 
turn again to this remark later in the analysis of the ®Li 
and ®Be EDMs. 



r(fm) 


FIG. 1: Folding of the CP-odd one pion exchange a — N 
potential. The CP-odd N — N interaction is also shown for 
comparison. The coupling constant was factored out. 


IV. RESULTS AND DISCUSSION 

We now give the result of the evaluation of the EDM of 
the ®Li nucleus in the cluster approximation. The result 
is shown in Table HI For the isovector CP-odd one-pion 
exchange contribution, we have found that the ®Li EDM 
enhances the CP violation by 2 times compared with that 
of the deuteron. The ®Li EDM is made of two compo¬ 
nents, namely the deuteron cluster polarization and the 
effect of the CP-odd a — TV interaction. The contribu¬ 
tion from the deuteron subsystem to the ®Li EDM takes 
a value close to the single deuteron EDM (about 43% of 
the total ®Li EDM). This can also be seen by inspect¬ 
ing the ®Li EDM induced by the isovector exchange of 
p and 77 mesons, which has no contribution to the CP- 
odd a — N potential. As the sensitivity to them takes 
values very close to those of the deuteron EDM, we may 
say that the deuteron cluster is relevant in the ®Li EDM. 
The CP-odd a — TV contribution to the ®Li EDM takes a 
comparable value as the contribution from the deuteron 
cluster. This additional effect makes the enhancement of 
the ®Li EDM. 

For the w meson exchange CP-odd potential, the en¬ 
hancement of the CP-odd effect in the ®Li is more accen¬ 
tuated. We must however be careful on it, since it may be 
an overestimation, as the short range part was averaged 
due to the folding. The calculated contributions of the 
isovector CP-odd p and r] exchange potential are reliable 
since their effects arise only from the deuteron cluster. 
We must also note that the isoscalar and isotensor CP- 
odd potentials do not contribute to the ®Li EDM, due to 
the inert a cluster and the isospin symmetry. It is how¬ 
ever to be noted that this is not the case if we consider 
the isospin breaking effect together with the excitation 
of the a cluster, since the spin and isospin structures be¬ 
come active [^. These topics are the subject of a more 
microscopic level study, and need to be investigated in 
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future works. 

We also give the result of the calculation of the ®Be 
EDM in the cluster approximation. It is shown in Table 
m The sensitivity of the ®Be EDM to the isovector CP- 
odd pion exchange nuclear force is comparable to that of 
the deuteron. Interestingly, the effect of the polarization 
due to the CP-odd a — N interaction is close to that for 
the ^Li EDM. 

For heavier meson ( 77 , p, w) exchange processes, the sen¬ 
sitivity of the ®Be on the CP-odd potential is much more 
important, but this is probably an overestimation, since 
all those effects arise from the dangerous folding of the 
heavy meson exchange CP-odd nuclear force. As for the 
®Li, there are also CP-odd interactions which may con¬ 
tribute with the excitation of the a cluster. A more care¬ 
ful inspection of them are therefore needed at a more 
microscopic level. 

We have also tested the ab initio calculation of the 
EDMs of the deuteron, ^He and using the Argonne 
7;18 interaction [s^ and the CP-odd hamiltonian of Eq. 

The results are shown in Table ID For the deuteron, 
the result is in good agreement with those of Refs. 
[13, 113 . We also point that the coefficients of the in¬ 
trinsic nucleon EDM contribution {<yp^n)^ii are smaller 
than one. In the literature, these values were often set 
to one, but this is only true if the deuteron total angular 
momentum does not receive any orbital angular momen¬ 
tum contributions. As the deuteron is known to have a 
d-wave component, the formula should be corrected. 

Let us now present the result for the ^He and nu¬ 
clear EDM. In our framework, the binding energies of 
the ^He and nuclei obtained with the Argonne tilS 
nuclear force are 6.93 MeV and 7.63 MeV, respectively. 
The difference from the experimental values (7.7 MeV for 
^He and 8.4 MeV for ^H) is understood as the lack of the 
three-body force contribution [d^l. For the EDMs of the 
^He and ^H, we have found values which are in agree¬ 
ment with the recent study using the chiral effective held 
theory [13, [13| • It is also in agreement with Ref. [13 for 
the CP-odd isoscalar and isovector nuclear forces. Our 
result differs from that of Ref. [13 by a factor of 1/2 for 
all CP-odd nuclear force contributions. 

Let us hnally see the prospects for the observation 
of the new physics BSM. If we model the new physics 
contribution by the exchange of new particles with mass 
Mtvp in the virtual state with 0(1) CP phase, the dimen¬ 
sional analysis gives the typical CP-odd N — N coupling 

as Gt, 9np1w^^ where Aqcd ~ 200 MeV, and gNP 
is the coupling between quarks and new particles. If the 
EDM of the ®Li nuclei can be measured at the level of 
O( 10 “^®)e cm, we thus can probe the new physics scale 
Mprp ^ PeV [with g^p = 0(0.1)]. This naive estima¬ 
tion works for models which generate isovector CP-odd 
4-quark interactions, such as the Left-right symmetric 


model ( 43 . 

For the supersymmetric model, the sensitivity of 
O(10“^®)e cm for the ®Li EDM can probe the CP phases 
of the g term 9^ and the trilinear supersymmetry break¬ 
ing coupling 9a at the level of 0(10“^) for the su¬ 
persymmetry breaking scale Msusy ~ TeV H. Here 
we have assumed the Peccei-Quinn symmetry [^, and 
tan/3 = 0(1). The sensitivity to 9^ may be increased 
with growing tan /3 @,[ 53 . This prospective sensitivity 
can thus unveil the CP violation in the high scale super- 
symmetry breaking scenario. 

We also give here the sensitivity to the class of mod¬ 
els which generate Barr-Zee type diagrams. This is 
the case for the Higgs doublet models JsJ, supersym¬ 
metric models with R-parity violation [3, l52ll. etc, and 
they contribute through the chromo-EDM [^. Using 
the simple formula of the quark chromo-EDM ~ 

111 , where ttiq is the inner loop quark mass, 

Yq and Yq are the couplings between the exchanged 
scalar and the quarks q and Q, respectively, we obtain the 
sensitivity on the scale of new physics Mjvp ^ 

PeV. If the coupling constants Yq and Yq are small, the 
sensitivity to the scale BSM is attenuated. 


V. SUMMARY 

In this paper we have calculated the EDM of ^H, ^He, 
^H, ®Li, and ®Be nuclei using the Gaussian expansion 
method. We have found that the ®Li enhances the CP- 
odd effect, due to the effect of CP-odd a — N interac¬ 
tion, in addition to the polarization contribution from 
the deuteron subsystem. With the experimental sensi¬ 
tivity of 0(10“^®)e cm for the ®Li EDM, we have strong 
chances to unveil the CP violation BSM at the PeV scale. 
We therefore strongly recommend the experimentalists to 
study and measure the EDM of the ®Li nucleus. 
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